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Abstract 

We theoretically analyze non-local effects in electron transport across three-terminal ballistic normal-superconducting- 
normal (NSN) structures with spin-active interfaces. Subgap electrons entering S-electrode from one N-metal may 
form Cooper pairs with their counterparts penetrating from another N-metal. This phenomenon of crossed Andreev 
reflection is highly sensitive to electron spins and yields a rich variety of properties of non-local conductance which 
we describe non-perturbatively at arbitrary interface transmissions, voltages and temperatures. Our results can be 
applied to hybrid structures with normal, ferromagnetic and half-metallic electrodes and can be directly tested in 
future experiments. 
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1. Introduction 

Andreev reflection (AR) [1] is the main mecha- 
nism of low energy electron transport between a nor- 
mal metal and a superconductor. This mechanism 
results in a number of interesting effects causing, 
e.g., a non-zero subgap conductance [2] of such hy- 
brid structures. In systems with one superconduct- 
ing (S) and several normal (N) terminals, e.g. NSN 
hybrid structures, electrons may suffer Andreev re- 
flection at each of NS interfaces. Provided the dis- 
tance between two NS interfaces L strongly exceeds 
the superconducting coherence length £, AR pro- 
cesses at these interfaces are independent. If, how- 
ever, the distance L becomes comparable with £, 
two additional non-local processes should be taken 
into account (see Fig. 1). Firstly, an electron with 
subgap energy can directly penetrate from one N- 
metal into another through a superconductor. Since 
the subgap density of states in the superconductor 
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vanishes, the probability of this process should be 
suppressed by the factor ~ exp(— £/£). Secondly, an 
electron penetrating into the superconductor from 
the first N-terminal may form a Cooper pair together 
with another electron from the second N-terminal. 
In this case a hole goes into the second N-metal 
and AR becomes a non-local effect. The probability 
of this process - usually called crossed Andreev re- 
flection (CAR) [314] - also decays as ~ exp(— L/£) 
and, in combination with direct electron transfer be- 
tween normal electrodes, determines non-local con- 
ductance in hybrid multi-terminal structures. This 
non-local conductance can be directly measured in 
experiments and it has recently become a subject of 
intensive investigations. 

Several experiments [5I6|7] allowed to clearly de- 
tect the non-local conductance in three-terminal 
NSN structures and demonstrated a rich variety of 
different features which unambiguous and detailed 
interpretation remains an important task for the fu- 
ture. At this point we note that in addition to CAR 
a number of other physical effects may consider- 
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Fig. 1. Two elementary processes contributing to non-local 
conductance of an NSN device: (1) direct electron transfer 
and (2) crossed Andreev reflection. 



ably influence the observations. Among such effects 
we mention, e.g, charge imbalance (relevant close 
to the superconducting critical temperature |5l7j ) 
and zero-bias anomalies in the Andreev conduc- 
tance due to both disorder-enhanced interference of 
electrons [819)10] and Coulomb effects |10|ll|12j . 

Theoretically CAR was analyzed within the per- 
turbation theory in the transmission of NS inter- 
faces in Refs. [13114) where it was demonstrated that 
in the lowest order in the interface barrier trans- 
mission and at T = CAR contribution to cross- 
terminal conductance is exactly canceled by that 
from elastic electron cotunneling (indicated as (1) in 
Fig. 1), i.e. the non-local conductance turns out to 
vanish in this limit. Recently a non-trivial interplay 
between normal reflection, tunneling, local AR and 
CAR in three-terminal ballistic NSN devices was 
non-perturbatively analyzed to all orders in the in- 
terface transmissions [15] . This analysis allowed to 
determine an explicit dependence of the non-local 
conductance both on the transmissions of NS inter- 
faces and on the length L, to set the maximum scale 
of the effect and to consider various important lim- 
its. The effect of disorder on CAR was recently stud- 
ied in Refs. [16] (perturbatively in tunneling) and 
[17j (non-perturbatively in tunneling, for a device 
with normal terminals attached to a superconduc- 
tor via an additional normal island). The interplay 
between CAR and Coulomb interaction effects was 
recently addressed in Refs. [18119] . 

It is also important to mention that both AR 
and CAR should be sensitive to magnetic proper- 
ties of normal electrodes because these processes 
essentially depend on spins of scattered electrons. 
First experiments on ferromagnet-superconductor- 
ferromagnet (FSF) structures [5] illustrated this 
point by demonstrating the dependence of non- 
local conductance on the polarization of ferromag- 




Fig. 2. Schematics of our NSN device. 



netic terminals. Hence, for better understanding of 
non-local effects in multitcrminal hybrid proxim- 
ity structures it is desirable to construct a theory 
of spin-resolved CAR. In the lowest order order in 
tunneling this task was accomplished in Ref. [13] . 
For FSF structures higher orders in the interface 
transmissions were considered in Refs. [20121] . 

In this paper we are going to generalize our qua- 
siclassical approach [15] and construct a theory of 
spin-resolved CAR to all orders in the interface 
transmissions. Instead of dealing directly with FSF 
devices we will consider NSN structures with spin- 
active interfaces. This model allows to distinguish 
spin-dependent contributions to the non-local con- 
ductance and to effectively mimic the situation of 
ferromagnetic and/or half-metallic electrodes. 

The structure of the paper is as follows. In Sec. 
2 we will introduce our model and discuss the qua- 
siclassical formalism supplemented by the bound- 
ary conditions for Green-Keldysh functions which 
account for electron scattering at spin-active inter- 
faces. Non-local electron transport in NSN struc- 
tures with such interfaces will be analyzed in Sec. 3. 
Our main conclusions will be briefly summarized in 
Sec. 4. Technical details related to boundary condi- 
tions will be outlined in Appendix A. 

2. The model and formalism 

Let us consider three-terminal NSN structure 
depicted in Fig. [2] We will assume that all three 
metallic electrodes are non-magnetic and ballistic, 
i.e. the electron elastic mean free path in each metal 
is larger than any other relevant size scale. In order 
to resolve spin-dependent effects we will assume 
that both NS interfaces are spin-active, i.e. we will 
distinguish "spin-up" and "spin-down" transmis- 
sions of the first (-Dij and D14) and the second (Z?2T 
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and D21) SN interface. All these four transmissions 
may take any value from zero to one. The effective 
cross-sections of the two interfaces will be denoted 
respectively as A\ and A 2 . The distance between 
these interfaces L as well as other geometric param- 
eters are assumed to be much larger than yMi^j, 
i.e. effectively both contacts are metallic constric- 
tions. In this case the voltage drops only across SN 
interfaces and not inside large metallic electrodes. 
Hence, nonequilibrium (e.g. charge imbalance) ef- 
fects related to the electric field penetration into 
the S-electrode can be neglected. In what follows 
we will also ignore Coulomb effects [10111112] . 

For convenience, we will set the electric potential 
of the S-electrode equal to zero, V = 0. In the pres- 
ence of bias voltages V\ and V 2 applied to two nor- 
mal electrodes (see Fig. ^ the currents I\ and I 2 
will flow through SNi and SN2 interfaces. These cur- 
rents can be evaluated with the aid of the quasiclassi- 
cal formalism of nonequilibrium Green-Eilcnberger- 
Keldysh functions g RAK [22] which we briefly spec- 
ify below. 

2.1. Quasiclassical equations 

In the ballistic limit the corresponding equations 
take the form 

ef 3 + eV{r, t) - A(r, t), g R ^ K {p F , e, r, t) 

+iv F Vg R ^ K [p Fl e 1 r 1 t) = 0, 

(1) 

where [a, b] = ab — ba, e is the quasiparticle energy, 
Pf = mvp is the electron Fermi momentum vec- 
tor and T3 is the Pauli matrix in Nambu space. The 
functions g R > A > K also obey the normalization con- 
ditions [g R f = (g A ) 2 = 1 and g R g K + g K g A = 0. 
Here and below the product of matrices is defined 
as time convolution. 

Green functions g R > A ' K and A are 4x4 matrices 
in Nambu and spin spaces. In Nambu space they can 
be parameterized as 



-R,A,K 



^R,A,K jR,A,K 
JR,A,K ~R,A,K 
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where g R > A > K , f R > A , K , f R , A ,K^ gR,A,K are 2 x 2 ma- 
trices in the spin space, A is the BCS order param- 
eter and ai are Pauli matrices. For simplicity we 
will only consider the case of spin-singlet isotropic 
pairing in the superconducting electrode. The cur- 



rent density is related to the Keldysh function g K 
according to the standard formula 

3{r,t) = - 



J de(v F Sp[T 3 g K {p F ,e,r,t)}) , 



( 3 ) 

where iVo = mp F /27r 2 is the density of state at the 
Fermi level and angular brackets (...) denote aver- 
aging over the Fermi momentum. 

2.2. Riccati parameterization 

The above matrix Green-Keldysh functions can 
be conveniently parameterized by four Riccati am- 
plitudes ry R - A , j R > A and two "distribution func- 
tions" x K , x K (here we follow the notations adopted 
in Ref. [23]): 



) K = 2N R 



X K _ ^K-A 



j R x K +x K ^ A x K -7 R x K ~/ A ) ' 

(4) 

where functions j R > A and j R ' A are Riccati ampli- 
tudes 



~R,A 



= ±N 
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-2j R < A 



2j R > A 

_ 1 _^R,A^R,A 
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N R ' A = 



and N R ' A are the following matrices 

(i_^ 7 «.^)-i 

'(6) 

With the aid of the above parameterization one can 
identically transform the quasiclassical equations 
([1]) into the following set of effectively decoupled 
equations for Riccati amplitudes and distribution 
functions [23j 

iv F V~f R > A + [e + eV(r, t)}~/ R ' A + j R ' A [e - eV(r, t)] 

= 7 i ^A*ia 27 fl < A - Aia 2l 
(7) 

iv F Vj R ' A - [e - eV(r, t)]j R - A - j R < A [e + eV(r, t)] 

= ^ A Ma 2 ^ A - A*ia 2 , 
(8) 

iv F \7x K + [e + eV(r, t)]x K - x K [e + eV(r, t)} 

-y R A*ia 2 x K - x K Aia 2 j A = 0, 

(9) 

iv F Vx K - [e - eV(r, t)]x K + x K [e - eV(r, t)} 
-j R Aio- 2 x K ~ x K A*ia 2 -y A = 0. 

(10) 
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Fig. 3. Riccati amplitudes for incoming and outgoing trajectories for an NSN structure with two barriers. The arrows define 
quasiparticle momentum directions. We also indicate relevant Riccati amplitudes and distribution functions parameterizing 
the Green-Keldysh function for the corresponding trajectory. 



Depending on the particular trajectory it is also 
convenient to introduce an additional "replica" of 
both Riccati amplitudes and distribution functions 
which - again following the notations adopted in 
Refs. |23|24j ) will be denoted by capital letters T 
and X . These "capital" Riccati amplitudes and dis- 
tribution functions obey the same (|T|)- (flT)]) with the 
replacement 7 — > T and x — > X . The distinction be- 
tween different Riccati amplitudes and distribution 
functions will be made explicit below. 

2.3. Boundary conditions 

Quasiclassical equations should be supplemented 
by appropriate boundary conditions at metallic in- 
terfaces. These conditions were initially formulated 
by Zaitsev [25] and later generalized to the case of 
spin-active (and specularly reflecting) interfaces in 
Ref. [28]. 

Before specifying these conditions we would like 
to note that the applicability of the above quasiclas- 
sical formalism with appropriate boundary condi- 
tions to hybrid structures with two (or more) bar- 
riers is, in general, a non-trivial issue |26|27] which 
requires a comment. Electrons scattered at different 
barriers may interfere and form bound states (res- 
onances) which in general cannot be correctly de- 
scribed within the quasiclassical formalism employ- 
ing Zaitsev boundary conditions. In our geometry, 
however, any relevant trajectory reaches each inter- 
face only once whereas the probability of multiple 
reflections at both interfaces is small in the param- 
eter A1A2 1 L 4 <C 1. Hence, resonances formed by 
multiply reflected electron waves can be neglected, 



and our formalism remains adequate for the prob- 
lem in question. 

In what follows we will make use of boundary con- 
ditions formulated directly in terms of Riccati am- 
plitudes and the distribution functions. In the case 
of spin-active interfaces these conditions are rather 
cumbersome [23] and therefore are relegated to Ap- 
pendix. 

In our three terminal geometry nonlocal conduc- 
tance arises only from trajectories that cross both 
interfaces. Consider first the contribution of trajec- 
tories illustrated in Fig. [3] where, for simplicity, we 
assume identical polarizations of both interfaces SNi 
and SN2. In this case the scattering matrix S for 
the first interface is defined in Eq. (|A.23|) while the 
5-matrix for the second interface can be obtained 
from Eq. (|A.23j) by means of the replacement 1 — > 
2. Accordingly, all Riccati amplitudes have zero di- 
agonal elements and all distribution functions have 
zero off-diagonal elements. 

Finally, one needs to specify the asymptotic 
boundary conditions far from NS interfaces. Deep 
in metallic electrodes we have 
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(14) 

while in the bulk of superconducting electrode we 
have 



it 



= — tanh 



2T 
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7^ = -a(e)ia 2 , 7i' = a*(e)ia 2 , 
[1-K £ )| 2 ]tanh^, 



2T' 

7 ^ = a(e)ia 2 , J§ = -a*[e)ia 2 
.x£ = [l-|a( £ )| 2 ]tanh^, 



(15) 
(16) 
(17) 
(18) 



where a(s) = — (e — \/ e 2 — A 2 )/A. 



2.4. Green functions and non-local conductance 

With the aid of the above equations and boundary 
conditions it is straightforward to evaluate the qua- 
siclassical Green-Keldysh functions for our three- 
terminal device along whole trajectory. For instance, 
from the boundary conditions on the second inter- 
face we find 



r|5 — y/ R 2 ^R 2 l 



J02 







e 



ia 2 a(e). (19) 



Then integrating equation ((7]) along trajectory con- 
necting both interfaces and using Eq. (TIT?)) as the 
initial condition we arrive at the Riccati amplitude 
at the first interface: 



7? 



T§ + (eT R , + Aia 2 )Q 
l-(e-AT§ia 2 )Q ' 



tanh [iQ,L/vf\ n 
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Similarly, integrating Eq. ^ , one finds 

1 - 



x. 
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i-r*r£ 
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We also note that the relation (7-^)+ = f A ' R 
makes it unnecessary (while redundant) to sepa- 
rately calculate the amplitudes 7j; and Tf, . 

Finally we arrive at the following general expres- 
sion for the Green-Keldysh function g^ at SNi in- 
terface UHl 



9? 



5n(Vi) + ^f 2 (^i) + ^f(^) 



(22) 



where the three different terms in the right-hand 
side correspond respectively to local contribution 
and two non-local corrections depending on V\ and 
V 2 . Since here we are only interested in the non- 
local conductance of our structure, it is sufficient to 
keep only the part of the Green-Keldysh function 
Sgf which depends on the voltage V 2 across the sec- 
ond interface. Combining the above results with the 
boundary conditions for the distribution functions 



and X^ at both interfaces respectively one im- 
mediately recovers Sg^ . E.g., for the outgoing mo- 
menta directions (l out ) one gets 
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5Xf = 2(1 - tanh 2 im/v F ) 
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where SX^ is the V^-dependent part of the X^ dis- 
tribution function, Zj = ^JFL^R^ exp(i#j), i = 1, 2 
Ri-\(D are the spin-sensitive reflection coefficients, 9i 
are the spin-mixing angles and 

2 tanhiLfi/up 



P(e,z 1 ,z 2 ) 



1 — z\z 2 a 



n 



x [e(l + ziz 2 a 2 ) + aA(z 1 + z 2 )] 



(24) 



Reversing momenta directions for all trajectories 
in Fig. [3] yields the second set of trajectories which 
also contributes to the non-local conductance. The 
corresponding expression for the V 2 -dependent part 
of the Green-Keldysh function bgf^ at the nor- 
mal metal side of the first interface is derived analo- 
gously to Eq. (|23|) . Again for the outgoing momenta 
directions we obtain 



<5gf (2) = 2(1 - tanh 2 im/v F ) 



l/P(e,Zi,z 2 ) 


Dx^RiiD 2 ^R 2 i 
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Finally we note that in the expression for function 
g K there also exist terms corresponding to incoming 
momenta directions (again on the normal metal side 
of the first interface). However, those terms do not 
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depend on the voltage V2 and, hence, do not con- 
tribute to the non-local conductance at all. Thus, 
Eqs. (|23|) and (|25|) already contain complete infor- 
mation which allows to evaluate the non-local cur- 
rent across the first interface. This task will be ac- 
complished below. 

3. Nonlocal conductance 

3.1. General results 

With the aid of general expressions for the Green- 
Keldysh functions obtained here it is possible to 
demonstrate that the total currents across the first 
and the second interface NS interfaces I\ and I2 have 
the form [T5l 



h=hl(V 1 )+I 1 2{V 2 ), 
h=I 2 l{V 1 )+I 22 {V 2 ). 



(26) 
(27) 



Here we are interested in the non-local part of the 
current 1x2^2) across the first interface. Substitut- 
ing $23j) and d25j) into © we obtain 
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de 
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where 
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JV i2 



Go 



D\]D 2 ] + D\iD 2 i 



(29) 



is the non-local conductance of our device in the 
normal state, 

87172A/1W2 



G = 



Pf7i(2) is normal to the first (second) interface 
component of the Fermi momentum for electrons 
propagating straight between the interfaces, N\^ 2 = 
p 2 F A\ y2 / Ait define the number of conducting chan- 
nels of the corresponding interface, R q = 2ir/e 2 is 
the quantum resistance unit. 

Eq. (|28|) is the central result of our paper. It fully 
determines the non-local spin-dependent current in 
our three-terminal ballistic NSN structure at arbi- 
trary interface transmissions, voltages and tempera- 
ture. It is also convenient to introduce the non-local 
differential conductance Gi 2 (V 2 ) = —dI\ 2 jdV 2 . In 



the limit T,V±. 2 <C A only subgap quasiparticles 
contribute to the current and the differential con- 
ductance becomes voltage-independent. We have 
1 12 = -G12V2, where 

G12 = G (l - tanh 2 LA/v F ) 

(git - guKggT - D21) + D^DnD^Dn 
P(0,z 1 ,z 2 ) 

(30) 

For spin- isotropic interfaces Eqs. (|3TJ[) and (1281) yield 
the results obtained previously in [T5] . In the lowest 
(first order) order in the transmissions of both in- 
terfaces Eq. (|3"0|) reduces to the result by Falci et al. 
[T3] provided we assume that the ratio Du 2 u /Dx( 2 )l 
coincides with that of spin-up and spin-down densi- 
ties of states in the first (second) ferromagnetic elec- 
trode. We also note that, provided at least one of 
the interfaces is spin-isotropic, the conductance (|30p 
is proportional to the product of all four transmis- 
sions Di^D\iD2]D2[, i.e. it differs from zero only 
due to processes involving scattering with both spin 
projections in both normal electrodes. 

As in the case of spin-isotropic interfaces the 
value G12 ((317)) gets strongly suppressed with in- 
creasing L, and this dependence on L is in general 
non-exponential reducing to G12 oc exp(— 2LA/vp) 
either in the limit of small transmissions or large 
L > v F /A. 

In the spin-degenerate case for a given L the non- 
local conductance reaches its maximum reflection- 
less barriers -D1.2 — 1- Interestingly, in this case for 
mall L <§; vf/A the conductance G12 identically 
coincides with its normal state value Gj\r 12 at any 
temperature and voltage. For L — > there is "no 
space" for CAR to develop on these trajectories and, 
hence, CAR contribution to G12 vanishes, whereas 
direct transfer of electrons between Ni and N2 re- 
mains unaffected by superconductivity in this limit. 

The situation changes provided at least one of the 
transmissions is smaller than one. In this case scat- 
tering at SN interfaces mixes up trajectories con- 
necting Ni and N 2 terminals with ones going deep 
into and coming from the superconductor. As a re- 
sult, CAR contribution to G12 does not vanish even 
in the limit L — > and G12 turns out to be smaller 
than Gat 12 . 

3.2. Polarized interfaces 

Let us now turn to the limit of highly polarized 
interfaces which is accounted for by either spin-up 
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Fig. 4. Zero temperature differential non-lo- 
cal conductance as a function of voltage for a 
half-metal/superconductor/half-metal structure. Different 
curves correspond to different interface transmissions and 
different orientations of the half-metal magnetizations. 

or spin-down transmission of each interface going 
to zero. In this limit our model describes an HSH 
structure, where H stands for fully spin-polarized 
half-metallic electrodes. In the case of parallel mag- 
netization of both half-metals (Di^ = D\, = 0, 
D 2 ] = D 2 , and D21 = 0) we obtain 



12(^2) = 
x(l 



Go 
4^ 



de 
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P(e,z 1 ,z 2 ) 

while in the case of antiparallel magnetization 
(D 1T = D u D n = 0, L» 2T = 0, and D 2l = D 2 ) we 
find 



tanh ■ 



x (1 - tanh 2 iLtt/v F ) 



— — tanh — 

2T 2T 

D x D2%a\; 



(32) 



P(-£, Zl, Z 2 ) 

i.e. the non-local conductances |3T|) and (|32|) have 
opposite signs. Some typical curves for the differen- 
tial non-local conductance are presented in Fig. [4] 
at sufficiently high interface transmissions and zero 
spin mixing angles. With decreasing interface trans- 
missions low voltage non-local conductance dimin- 
ishes and sharp peaks at voltages eV 2 — ±A appear. 

We will now turn to asymmetric HSN het- 
erostructures with one half-metallic and one spin- 
isotropic electrode. In this case the non-local cur- 
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Fig. 5. Zero temperature differential non-lo- 
cal conductance as a function of voltage for a 
half-metal/superconductor /metal structure. Different curves 
correspond to different values of the spin-mixing angle d\ . 

rent does not depend on the magnetization di- 
rection of the half-metal. Below we will distin- 
guish NSH and HSN configurations. In the case of 
metal/superconductor/half- metal (NSH) structures 
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We note that in this case the spectral cur- 
rent is identically zero at all subgap energies 
|e| < A. The sign of the non-local current and 
differential conductance remain the same as in 
the normal state. In the reversed situation of a 
half-metal/superconductor/metal (HSN) structure 



(D n =Di,D n 
results yield 



0, D 2 f = D21 = D 2 ) our general 
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P(e, zi,z 2 ) P{—e,zi,z 2 )J 
In contrast to Eq. (|33|) . here the spectral current 



does not vanish at subgap energies provided the 
spin-mixing angle 9i differs from zero. The sign of 
the non-local conductance in this geometry depends 
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Fig. 6. Zero temperature differential non-local conductance 
as a function of voltage for a general case of an NSN structure 
with non-zero transmissions , Di\ and D21 and 

non-zero spin-mixing angles 8± and #2 ■ 

on both applied voltage V% and spin-mixing angle 9\ . 
The zero temperature non-local conductance evalu- 
ated with the aid of Eq. at sufficiently high in- 
terface transmissions and different spin-mixing an- 
gles is depicted in Fig. 

In general, the voltage dependence of the non- 
local current is very sensitive to the interface trans- 
missions and - in particular - to values of spin- 
mixing angles 9\^- Typical curves illustrating the 
voltage dependence of the differential non-local con- 
ductance are shown in Fig.[6]at different spin-mixing 
angles. Positions of the peaks correspond to energies 
of the quasibound states in NSN structures. 

4. Conclusions 

In this paper we have developed a non-perturbative 
theory of spin-resolved non-local electron transport 
in ballistic NSN three-terminal structures with spin- 
active interfaces. Our theory applies at arbitrary in- 
terface transmissions and allows to fully describe a 
non-trivial interplay between spin-sensitive normal 
scattering, local and non-local Andreev reflection at 
SN interfaces. We have evaluated the non-local con- 
ductance of our NSN device at arbitrary voltages 
and temperatures and observed a number of inter- 
esting properties of such structures with spin-active 
interfaces. Our results can be applied to various 
NSN hybrid structures, including systems with fcr- 




Fig. A.l. Riccati amplitudes for incoming and outgoing tra- 
jectories from the both sides of the interface. 

romagnetic and half-metallic electrodes, and can be 
directly tested in future experiments. 

Appendix A. Boundary conditions for 
spin-active interfaces 

Let us explicitly specify the relations between Ric- 
cati amplitudes and distribution functions for in- 
coming and outgoing trajectories at the first inter- 
face, see Fig. IA.ll The boundary conditions for Tf , 
rf and Xi can be written in the form [24] 

rf = T%$&+i%fi§&. t (A.i) 

f^=£ n 7iV, + £ii'7^ (A.2) 
Xi — r u x x f lr + t\iXyt\ r — auXi> o, lr . (A. 3) 

Here the transmission (t), reflection (r), and branch- 
conversion (a) amplitudes are defined as follows: 

(A.4) 
(A.5) 

rt = i)" 1 - SMPvv)- 1 }- 1 , 

(A.6) 

t? r = -iPvvr^iPtn)- 1 ~ Sii'OtfiO" 1 ]" 1 . 

(A.7) 

< = (T?S n - Sin?)0&,)-\ (A.8) 

< = (Pnr'iStJf - 7i A S+), (A.9) 

where 

oR_a+_ Ra+~R nR _ o _^Ro..^R 
Pij ~~ ij ij —ijli ' Pij ~ lj D 3 l li > 

(A.10) 

oA _ q /In ~A a A _ c+ -i c+ A 

Pij — °ij li i Pij ~ —ji 7« °jilj • 

(All) 

Similarly, the boundary conditions for Tf, Tf, and 
Xf are: 



8 



R z,R < 



1 1 — D n7i ' ir ' D x'X n'H 



yK _ ~R~K A , fn~i\ f ^ 
^1 — 'lr tl ll J 'l' l lr 



a l/ a 'l' a lr' 



(A.12) 
(A.13) 
(A.14) 



where 



zR 



(A.15) 
(A.16) 



(A.17) 

^r = -(^/iO- 1 tsiiG^i)- 1 -^i(^v)- 1 ]- 1 , 

(A.18) 

a5 = (rf5 n -5 ll7l ii )(^)- 1 ; (A.19) 

(A.20) 



Boundary conditions for T lt ' , f v ' , Xy and Xy 
can be obtained from the above equations by means 
of the replacement 1 «-> V . 

Matrices Six, Sw , Svi, and Svx' are the compo- 
nents of the iS-matrix describing electron scattering 
at the first interface: 



S 




SS^ 



(A.21) 



For specularly reflecting interfaces with inversion 
symmetry elements of the S- matrix have the follow- 
ing form 



5*11 — <5i'i' 



Sii' — Sn 




iSi/2 



(A.22) 



(A.23) 

where and are interface transmission for 
spin-up and spin-down electron. The spin-mixing 
angle 9\ accounts for the difference between scat- 
tering phases for processes with opposite spin direc- 
tions. The condition = Sij holds provided there 
exists reflection symmetry in the plane perpendicu- 
lar to the interface. 
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